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DOCUMENTATION 


Computation of the Total arid the B-avetage Bhattacharya Distance ; 

(Univac 1108, Univ. of Houston). » 

This program consists of 3 subroutines to be exfecuted in the following 
sequence: 

.(1) Subroutine BHATT 

(2) Subroutine BHATBl 

(3) Subroutine BHATB2 

1. SUBROUTINE BHATT 

ABSTRACT ' ^ ^ ^ ^ ^ ^ ^ 

This subroutine calculates the total Bhattacharyya Distance, BDIST, using 
all N channels. The output of this program, BDIST, will be used in comparing 
the difference " BDIST where H^ is the B-average Bhattacharyya 

Distance computed in the subroutines BHATBl, BHATB2. 

User's Information: ^ ' 

(Double Precision Version Only) . 

In order to use this subroutine the following FORTRAN calling sequence 
must be given : 

CALL BHATT (COVAR, XMEAN,M,N, BDIST) ■ 

where : 

COVAR(input) is a real 3”dlmenslonal array (MxNxN) and contains 

the liN^N class covariance matrices (positive de- 
finite symmetric) used as input. 



3MEAN(input) Is a real 2-dimensional array (M<N.) and contains 

the' M -dimensional class mean vectors. 

M(input) is the no. of classes under consideration i.e. the 

no. of covariance matrices and mean vectors. 

N(input) .is the dimension of the covariance matrices and the 

' mean vectors. 

BDIST (output) is the value of the total Bhattacharyya Distance com- 

puted by subroutine BHATT. 

SUBROUTINES USED ; 

Subroutine BHATT in turn calls the following subroutines 

1. Subroutine MATMUL. This subroutine computes the product of 2 

^ matrices'. It calls subroutines SUPSUM and ORDER. 

2. Subroutine GHLSKY. This subroutine computes the inverse of a 

positive definite symmetric matrix. 

3. Subroutine DET. This subroutine computes the determinant of a 
positive definite symmetric matrix. 


NOTE: (1). The format statements for input, output are dependent upon the 

dimensions of the input data and corresponding adjustments have to be made to 
formats when different sets of data are run. 

-r(2). The variables declared in the DIMENSION statements have to similarly 
correspond to the dimensions of the input data. 


ALGORITHM : . : 

Subroutine BHATT computes the value of the total Bhattacharyya Distance 
using the covariance matrices and mean Vectors as Inputs. ^ 
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The total Bhattacharyya Distance, BDIST, is computed by the formula 

\ i 

1 m 

BDIST = i ^ H(i.1) 

® i=l 

where H(i,j), the interclass Bhattacharyya Distance. between classes i and 
j is given by 


H(l,j) - exp[-.-^ 6^^ +• Zj) An 


where 6..= u, - u, and u, is the mean vector corresponding to class i 

IJ X j 1 

and is the covariance matrix corresponding to class i. 

2. SUBROUTINE BHATBl ; 

ABSTRACT 

This subroutine attempts to calculate the minimum B-average Bhattacharyya 
Distance using 1 Householder transformation to construct the B-matrix. 


USER’S INFORMATION ; 

(Double Precision Version Only) 

In order to use this subroutine the following FORTRAN calling sequence must 
be given: 

CALL BHATBl (COVAR, XMEAN, M,N, K, ITE, ALPHA) 


where 

■■ ■ ■ ■' , -r' 

COVAR(lnput) is areal 3-dlmensional array (MN^N) containing 

i 

the M Nxn covariance matrices. 
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XMEAN (Input) 

M(lnput) 

N (Input) 

K(input) 

ITE(input) 

ALPHA(input) 


is a real 2-dlmenslonal array (MxN) and contains 
the M N-dimeosional mean vectors used as input, 
is the number' of classes' under . consideration (i.e. 
the no. of covariance matrices and mean vectors) . 
is the dimension of the covariance matrices and the 
mean vectors. 

is the number of rows desired in the transformation 

matrix B (which is Kxti) 

is 1 + (the no', of iterations required) 

is a varying paraaieter in the iteration formula. 


OUTPUT OF SUBROUTINE BHATBl 

This subroutine has the following output: 

1. The transformation matrix B (which has dimension Kx*^ corresponding 
to a particular value of the Householder generator F.* 

2. The value of the B-average interclass Bhattacharyya Distance 

Hg(i* j) f i. ~ If . . . f 5 j “ i+lf . . . ftt 

3. The N-dimensional F-vector which is the generator of the House- 

X 

holder transformation H = I-2FF used in constructing the B-matrix 
B = (Ij^|z)H. 

4. The value of the B-average Bhattacharyya Distance, Hg corresponding 
to the matrix B. 

3Hg 

5. The partial derivative vector which contains, the partial 

derivatives of Hg with respect to the vector F. 

*See 'ALGORITHM* 
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Subroutines Used 

The following subroutines axfe in turn called, by ..sub routine BHATBl: 

1. Subroutine MATMUL - calls SUPSUM and ORDER. 

2. Subroutine CHLSKY. 

3. Subroutine DET. 

ALGORITHM 

Subroutine BHATBl attempts to compute the minimum B-average Bhattacharyya 
Distance using one Householder transformation to compute the B-matrix. The 
B-average Bhattacharyya Distance is given by the formula 


_ m-1 m 

^ S 21 H (i,j) 

^ “ i=l j&l ® 


where 


Hg(i,j) 


= expt- - I An(|Z^ + 


1/2 1 ^ ,1/2' 


■..A ' ^ 

where ~ ^i ~ ® ® matrix of rank K 

of the form B = (iJz)H where H = I-2FF^, | (f| | 1. An. initial guess for 

X 1 I T 

F is taken to be F^ = and the corresponding matrix 


B = 


■■ ■ X 

= (I Z) (I-2F F ) is computed. The corfesponding value of 
K-J O O 


. m-1 m 

Hr " i 21 2H HrUJ) 
® “ $1 j^I+i ^ 


is also computed. 
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The steepest descent Iterator Is then applied to alter the value of F 

3H 

P 

where a is a varying parameter and is one of the Inputs to the program. 

is the partial derivative vector (derived analytically). The value of 

‘’"p 

' \l M 

is then normalized SO that | = 1. The B-matrix is recomputed with 

the new value of F. The corresponding value of Hg is computed. This procedure 

is repeated (ITE - 1) number of times (8 seems to be a good value for ITE) . 

Two points should be. noted: 

3H 

(1) . Whether ^0. 

(2) . Whether ~ ^ “ BDIST (the total Bhattacharyya Distance) is 

sufficiently small. 

The values of a and ITE (which are both inputs to this subroutine) . 
should be altered accordingly in order to achieve the above 2 objectives. 

The value of, F at which the mlnimtan value of H„ occurs is saved. Call 

D 

it FI. , 

3. Subroutine BHATB2 

Tbis subroutine attempts to compute the minimum B-average Bhattacharyya 
Distance using 2 Householder transformations. 

USER’S INFORMATION ! 

(Double Precision Version) 

(1) In order to use this subioutine the' following FORTRAN calling 
sequence must be given: 
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t 

CALL BHATB2(C0VAR, XMEAN, M, N, K, ITE, ALPHA) 

where 


COVAR*. XMEAN, M,N,K,ITB, ALPHA 


have the same meanings as in SUBROUTINE BHATBl. 

(2) This subroutine reads in the value of FI computed in the previous 
program (subroutine BHATBl), The data cards for FI should have 
the format 5F16.8 (e.g. if FI is 12-dimensional then FI is 

^ punched on 3 data cards; the first 2 cards contain 5 components 

of FI and the last card contains 2 components of FI). 

These data cards for FI are placed following the data cards for the 
covariance matrices and the mean vectors. 

(3) The value of FI that is read in is then used to compute the 

- T 

Householder transformation H^= I - 2F1F1 . The covariance matrices 
and the mean vectors i = 1, .,.,m are transformed into 
Hf^i^l H^y^. 

The number of Householder transforpiations by which the covariance matrices 
and the mean vectors have to be transformed is denoted by the variable 

IJ. 

For Subroutine BHATB2 we require one Householder transformation to obtain 
H^Z^Hj^ and H^y^. 



8 


The FORTRAN statements "IJ «= 1" appears after the comment; 

”C— — — IJ Eq. No. of Householder Transformations Required-^—". 

OUTPUT OF SUBROUTINE BHATB2 . ' ' 

1. The vector FI. which is the generator of the Householder transfor- 
mation H^ = I - 2E1F1^. 

2. Same as subroutine BHATBl* 

ALGORITHM ; 

Here each is replaced by and each is replaced by 

The B matrix is then taken to be B ~ (I^^j Z) (I-2FF^) , F - 1. An initial 
T 1 1 * 

guess for F, * * * * niade and the same procedure as in subroutine 

BHATBl is applied. The value of F = F2 at which the minimum value of 
occurs is saved. 

USING MORE THAN 2 HOUSEHOLDER TRANSFORMATIONS TO CONSTRUCT THE B-MATRIX; 

If more than 2 Householder transformations are required to compute the 

transformation matrix B i. e. if ~ H_ - BDIST is not small enough, then 

n B 

subroutine BHATB2 can be modified in the following way. For the B-matrix^^^ ^ 
requiring 3 Householder transforaations do the following: 

, (1) Place the data cards containing the vector F2 (computed in the 
previous program) following the data cards containing FI. 

(2) The statement following the comment "C. .. Ij Eq. NO. OF HOUSE- 
HOLDER TRANSFORMATIONS REQUIRED ..." should be "IJ =2" 

For j S 4 Householder transformations required in computing the B-matrix* 
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(1) the data cards for FI,.,., F(J-l) should be placed after the data 
cards for the covariance laatrices and mean vectors; 

(2) the statement "IJ » 2” should be changed to "IJ = (J-1)". 
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